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(Dated: February 7, 2008) 

The branching fraction of the r~ — > 3/i~2/i + ^ T decay (h = ir,K) is measured with the BABAR 
detector to be (8.56±0.05±0.42) x 10" 4 , where the first error is statistical and the second systematic. 
The observed structure of this decay is significantly different from the phase space prediction, with 
the p resonance playing a strong role. The decay t~ — > /i(1285)7r — iv, with the /i(1285) meson 
decaying to four charged pions, is observed and the branching fraction is measured to be (3.9 ± 0.7 ± 
0.5) x 10" 4 . 

PACS numbers: 13.35.Dx, 14.60.Fg 



The high-statistics sample of r pair events collected by 
the BABAR Collaboration allows detailed studies of rare 
decays of the r lepton. This Letter presents a measure- 
ment of the t~ — > Zhr2h + v r decay Jh = tv,K) from a 
sample of over 34,000 such decays The large data 
set allows a first look into the decay mechanism and the 
search for resonant structure of the t~ — > 3h~2h + v T de- 
cay mode. The best previous measurement of the r~ — > 
3hT2h + v r branching fraction is (7.7 ± 0.5 ± 0.9) x 10~ 4 , 
based on 295 events by the CLEO experiment |2|. 

Tau decays to one and three charged hadrons have been 
used to test the Standard Model, measure the masses of 
the t~ and v T , study the properties of low-mass reso- 
nances, test CP violation in the lepton sector, and search 
for new physics. Moreover, the semi-leptonic decays of 
the t lepton are ideal for studying strong interaction ef- 
fects (for example, see Ref. 0) as they probe the matrix 
element of the left-handed current between the vacuum 
and the hadronic state Q. The results presented in this 
Letter suggest that further studies will be possible with 
t~ — > ?>h~2h + v T decays. 

This analysis is based on data recorded by the BABAR 
detector at the PEP-II asymmetric-energy e + e~ storage 
ring operated at the Stanford Linear Accelerator Center. 
The data sample consists of 232.1 fb _1 recorded at center- 
of-mass energies (y/s) of 10.58 GeV and 10.54 GeV be- 
tween 1999 and 2004. With a luminosity-weighted cross 
section for r-pair production of <7 TT = (0.89 ± 0.02) nb 
[f| , this data sample contains approximately 400 million 
t decays. Monte Carlo simulation is used to evaluate 
the background contamination and selection efficiency. 
The r pair production is simulated with the KK2f Monte 
Carlo event generator [j| and the r decays modeled with 
Tauola [6j according to measured rates [7j ■ 

The BABAR detector is described in detail in Ref. 
Charged particle momenta are measured with a five-layer 
double-sided silicon vertex tracker and a 40-layer drift 
chamber inside a 1.5-T superconducting solenoidal mag- 
net. A calorimeter consisting of CsI(Tl) crystals is used 
to measure electromagnetic-shower energies, and an in- 
strumented magnetic flux return (IFR) is used to identify 
muons. 

Since r pairs are produced back to back in the e + e~ 
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center-of-mass frame, the event is divided into two hemi- 
spheres in the center-of-mass frame based on the plane 
perpendicular to the thrust axis from the tracks in the 
event. Each hemisphere is assumed to contain the decay 
products of a single r lepton. The analysis procedure se- 
lects events with one track in one hemisphere (tag hemi- 
sphere) and five tracks in the other hemisphere (signal 
hemisphere). All tracks are taken as pions unless iden- 
tified as an electron or muon. The total event charge is 
required to be zero. 

Charged particles are required to have momentum 
greater than O.lGeV/c in the plane transverse to the 
beam axis. The distance of the point of closest approach 
of the track to the beam axis must be less than 1.5 cm 
(Axy)- In addition, the z coordinate (along the beam 
axis) of the point of closest approach of the track must 
be within 10 cm of the z coordinate of the production 
point. 

The background from non-r sources (in particular, 
Bhabha scattering and two-photon production) is re- 
duced by requiring the magnitude of the thrust (T) of the 
event to be between 0.92 and 0.99. The ratio PT/-E m issing 
is also used to reduce the background from two-photon 
production, which tends to have low pt and high i? m i ss ing- 
The pt is the transverse component of the vector sum of 
the momenta of all the charged particles in the event and 
-^missing is the missing energy in the event. Events are 
retained if they satisfy the following criteria: 

{PT/E missing > 0.3 and 0.92 < T < 0.93) or 
(PT /^missing > 0.2 and 0.93 < T < 0.94) or 
(PT/Emissing > 0.1 and 0.94 < T < 0.95). 

There is no requirement on Px/^missing if the thrust is 
between 0.95 and 0.99. 

Furthermore, reduction of the non-r background is 
made by requiring that the track in the tag hemisphere 
be identified as an electron or a muon and that the mo- 
mentum of the track in the center-of-mass frame be less 
than 4 GeV/ c. Electrons are identified with the use of the 
ratio of energy measured by the calorimeter to track mo- 
mentum (E/p), the ionization loss in the tracking system 
(dE/dx), and the shape of the shower in the calorime- 
ter. Muons are identified by hits in the IFR and en- 
ergy deposits in the calorimeter consistent with the min- 
imum energy hypothesis. Residual background from mul- 
tihadronic events is reduced by requiring that there be at 
most one electromagnetic calorimeter cluster in the tag 
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hemisphere with energy above 0.05 GeV. Further, the 
total neutral energy in the tag hemisphere must be less 
than 1 GeV. 

Additional criteria are applied to the five track system 
in the signal hemisphere to reduce background from pho- 
ton conversions. The event is rejected if any of the tracks 
is identified as an electron or if any pair of oppositely 
charged tracks is consistent with originating from a pho- 
ton conversion. The invariant mass of the five charged 
particles is required to be less than 1.8GeV/c 2 . All in- 
variant masses shown are calculated assuming that the 
particles are pions. 

It is also required that there be no tt° candidates in 
the signal hemisphere. A ir° candidate consists of two 
clusters in the electromagnetic calorimeter that are not 
associated with any track. Each cluster is required to 
have an energy of at least 0.050 GeV and the two clus- 
ters have a combined invariant mass between 0.115 and 
0.150 GeV/ c 2 . In addition, any remaining clusters with 
energy greater than 0.5 GeV that are not associated to a 
track are considered a tt° candidate. 

A total of 20920 and 13929 events are selected when an 
electron or muon, respectively, are identified in the tag 
hemisphere. 

The selection efficiency is defined as the number of 
events with ar-> "ih~2h + v T decay in signal hemisphere 
and a tau leptonic decay in the tag hemisphere divided 
by the number of t pair events with a t~ — > 3h~2h + Vj-. 
The branching fraction of the r leptonic decay mode 
is incorporated into the selection efficiency. The efficien- 
cies are (4.71 ±0.05)% and (3.03 ±0.04)% in the electron 
and muon samples, respectively. The efficiencies are ob- 
tained from the Monte Carlo simulation and the quoted 
uncertainty is the Monte Carlo statistical error. 

The background in the selected sample comes from 
other t decays and multihadronic events. The back- 
ground percentages in the electron and the muon tag 
samples estimated from the Monte Carlo simulation are 
(20.6 ± 2.0)% and (21.7 ± 2.1)%, respectively. The er- 
rors are the combined statistical and systematic uncer- 
tainties. The sources of background in the electron tag 
sample can be broken down into the following categories: 
t~ — * 3h~2h + n°v T decays (7.2%), r decays with one 
or three tracks and at least one n° (6.3%), r decays 
with a Kg (4.9%), multihadronic events (1.8%, primar- 
ily cc events) and a residual amount from other t de- 
cays (0.5%). Background from Bhabha scattering and 
two-photon production is negligible. The relative uncer- 
tainties range between 15% and 20% for each background 
and reflect the statistical precision of the data and Monte 
Carlo samples used to evaluate the backgrounds. The 
backgrounds in the muon tag sample are very similar. 

In order to validate our Monte Carlo simulation for 
the background contamination we use experimental data 
samples where the particular background is enhanced. 
The uncertainty on the t~ — * 3h~2h + ir°is T background 
is estimated to be 20% by comparing the number of ir° 
mesons reconstructed in five charged track sample in the 



data and Monte Carlo simulation. The background from 
t~ — > h~(> lTr°)iy T and r~ — > h~h~h + (> 1tt°)v t arises 
when one or both of the photons from the decay of a 
ir° converts to an e + e~ pair or from a ir° — > e + e~7 de- 
cay. The uncertainty on this background is estimated to 
be 15% from the number of conversions and number of 
tracks identified as electrons. 

Background can also arise from t~ — ► ir~ K®K®v T and 
t~ — ► hr h~ h + K2,v T decays, both of which have been ob- 
served by other experiments \jy. The background from 
these decays is determined by fitting the mass distribu- 
tion of 7r + 7r~ pairs to obtain an estimate of the number of 
Kg mesons. The background estimation uses the Monte 
Carlo prediction for the t~ — * 7r~ K^K^v t decays modes. 
The t~ — * h~h~h + KgV r decay mode is not simulated 
and the background is assumed to be the excess of Kg 
mesons in the data over the Monte Carlo prediction. The 
uncertainty in the background from r decays with Kg 
mesons was found to be approximately 20% and includes 
contributions from the statistical uncertainties of the fits 
to the mass distribution of 7r + 7r~ pairs and the branch- 
ing ratios of the background decay modes. In addition, 
checks were made to ensure that the Kg background was 
from t decays and not multihadronic events. 

The background from multihadronic events was esti- 
mated from the number of events for which the recon- 
structed mass of the five tracks is above the r mass, and 
also from the number of events with more than one elec- 
tromagnetic cluster in the tag hemisphere. The uncer- 
tainty in the multihadronic background is estimated to 
be 20%. 

The branching fraction is defined as B — N sc \(l — 
/bkgd)/(2iVe) where N se \ is the number of selected events, 
N is the number of tau pair events determined from 
the cross section and luminosity, /bkgd is the fraction of 
background, and e is the efficiency for selecting r~ — > 
2ih~2h + v T and lepton events. 

The branching fraction of the r~ — > 3h~2h + v T decay 
is found to be (8.53 ± 0.06 ± 0.42) x 10" 4 and (8.73 ± 
0.07 ± 0.48) x 10~ 4 for the data selected by the electron 
and muon tags, respectively. The first uncertainty is the 
statistical error and the second systematic. The average 
branching fraction is (8.56±0.05±0.42) x 10" 4 where the 
correlation between the systematic errors in the electron 
and muon tag results is taken into account. Our value 
of the branching fraction is in good agreement with the 
Particle Data Group fit value of (8.2 ± 0.6) x 10~ 4 0. 

The systematic error includes contributions from the 
efficiency for reconstructing the six tracks in the event 
(3.1%), the background in the sample (2.4%), the lu- 
minosity and t + t~ cross section (2.3%), the tt° finding 
algorithm (2.0%), and the lepton identification in the tag 
hemisphere (1.0% for electrons and 2.5% for muons). 

The error on the efficiency for reconstructing a track is 
estimated to be 1.2% for tracks with pt < 0.3 GeV/c and 
0.5% for tracks with px > 0.3 GeV/c. The errors were ob- 
tained from comparison of efficiencies of the standalone 
track reconstruction in the silicon vertex tracker and the 
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drift chamber, and confirmed by an independent analysis 
of t decays into three charged particles and a neutrino. 
Variation of selection cuts such as the minimum trans- 
verse momentum of the track, the number of tracks with 
hits in the silicon vertex tracker, and the sum of the dxy 
of the five tracks resulted in a negligible change in the 
branching fraction. 

Variation of the selection criteria produced consistent 
results for the branching fraction. In addition, the selec- 
tion efficiency was found to have no dependence on the 
reconstructed mass of the five tracks. 

In Fig. ^ the distribution of the invariant mass of the 
five charged particles in the signal hemisphere is pre- 
sented. The discrepancy between Tauola, which uses a 
phase space distribution for t~ — * 3n~ 2n + v T and the 
data is believed to be due to resonant contributions in the 
t~ — > 37r~27r + ^ r decay mode. There are three allowed 
isospin states for this decay mode (see Ref. Q) and two 
of these isospin states have particles with quantum num- 
bers of the p meson. Fig.|21shows the mass of h + h~~ pah- 
combinations where the shoulder at 0.77 GeV/c 2 suggests 
a strong contribution from the p resonance. 

No attempt was made to extract the fraction of p 
mesons as no model for resonant structure of the r~ — > 
Zh~2h + v T decay exists. Such a model would need to in- 
clude the three allowed isospin states and the admixture 
of the isospin states could be extracted from this data 
sample as it was done for r~ — ► h~h~h + v T [To|. 
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FIG. 1: Invariant mass of the five charged particles in the 
signal hemisphere after all other selection criteria (except the 
mass requirement) are applied. The points are the data and 
the histogram is the Monte Carlo simulation for both the elec- 
tron and muon tag samples. The unshaded and two shaded 
histograms are the signal, tau and multihadronic background 
events, respectively. The arrow indicates the selection re- 
quirement applied to the samples. The Monte Carlo sample 
is normalized to the luminosity of the data sample. 

The data sample can also be used to study the t~ — > 
/i(I285)7r~i/ T decay, where the /i(1285) decays into a 
2ir~2ir + final state. In Fig. [31 the invariant mass of the 
2h + 2h~ particle system is plotted for data. The fit to the 
data uses a second-order polynomial distribution for the 
background and a Breit-Wigner for the peak region. The 
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FIG. 2: Reconstructed mass of h + h~ pairs in the five tracks 
in the signal hemisphere. The data are shown as points with 
error bars. The unshaded and shaded histograms are the sig- 
nal and background predicted by the Monte Carlo simulation. 
The peak at 0.5 GeV/c 2 is due to K® mesons that are not re- 
jected by the selection. There are six entries per event. 



Breit-Wigner is convoluted with a Gaussian distribution 
with a standard deviation corresponding to the expected 
mass resolution. The background distribution was deter- 
mined by fitting the region between 1.1 and I.4GeV/c 2 
excluding the /i(I285) peak (1.25-1.31 GeV/c 2 ). 

A total of 1369 ± 232 r~ -> /i(1285)7r _ i/ T decays 
are obtained from the fit. The fraction of t~ — » 
/i(1285)7r _ ^ r decays found in the r~ — > 2>h~2h + v T sam- 
ple is measured to be (0.050 ± 0.008 ± 0.005) and the 
branching fraction of the t~ — ► /i(1285)7r - z/ T decay is 
calculated to be (3.9 ± 0.7 ± 0.5) x 10~ 4 . The branching 
fraction for the /i(1285)— > 2ir~2ir + decay used to calcu- 
late the t~ — ► ,/i(1285)7r~iA r branching fraction is taken 
from the Particle Data Group 0. The first errors are 
the statistical uncertainties obtained from the fit and the 
second errors are the systematic uncertainties. The sys- 
tematic uncertainties include a contribution from the fit 
(10%) estimated by studying the results of fits using dif- 
ferent mass bins, background functions and detector res- 
olutions. The systematic error on the branching fraction 
also includes the uncertainty on the branching fractions 
of the r--> Zh-2h+v T (5%) and the /i (1285) -> 27T-27T+ 
decay modes (6%). 

Checks confirmed that the /i(1285) signal did not arise 
from multihadronic events. This was done by relaxing 
the selection criteria in a way which increased the mul- 
tihadronic background and confirming that the /i(1285) 
signal did not increase. In addition, the observation of 
the t~ — > /i(1285)7r _ i/ r decay was confirmed by looking 
at a data sample with a hadron tag. 

Our value of the t~ — > /i(1285)7r - i/ r branching frac- 
tion is in agreement with the result obtained by the 
CLEO Collaboration, (5.8 ± 2.3) x 10~ 4 , obtained using 
the /i(1285)— *• 777T7r decay mode ^l|. It is also consistent 
with a theoretical prediction of 2.91 x 10~ 4 [T^ ]. 

In summary, the BABAR Collaboration has measured 
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FIG. 3: Reconstructed mass of the 2h + 2h~ combinations in 
the signal hemisphere. The solid line is a fit to the data using 
a second-order polynomial distribution (dashed-line) for the 
background and a Breit-Wigner convoluted by a Gaussian for 
the peak region. The data are shown as points with error 
bars. There are three entries per event. 

the t" — > Zh^2h + v T branching fraction, B(t~ — > 
ih-2h+v T ) = (8.56 ± 0.05 ± 0.42) x 10~ 4 . The mass 
of the five charged hadron system is not well described 
by a phase space model. The invariant mass distribution 
of h + h~ pairs shows that the p meson is produced in the 
t~ — * 3h~2h + v T decay. The decay r~ — > /i(1285)7r _ ^ T 



is confirmed in the /i(1285)— > 27r 27t + channel and the 
branching fraction measured is B(t~ — > /i (1285)7r - 2/ T ) = 
(3.9 ±0.7 ±0.5) x 10~ 4 . 
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